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Abstract 

A simple model for the nucleation and growth of single wall carbon nanotubes from a graphene 
sheet at the surface of a metallic catalyst saturated in carbon is developed. It enables to predict 
the geometry and energy of tube embryos of all possible chiralities, as well as the way that they 
can grow. It is shown that armchair-like chiralities are energy preferred for geometrical reasons. 
This result is discussed and compared to experimental literature. 

1 Introduction 

Since their discovery twenty years ago, carbon nanotubes (CNT) have attracted a very strong interest 
among the scientific community. They are indeed fascinating objects on an academic point of view, and 
they are foreseen for a large panel of future applications. Among these potential applications, some 
rely strongly on the helicity of CNTs, mostly concerning electrical transport properties, but present 
elaboration techniques do not enable to control this helicity: this is a challenge for the future. It is 
highly desirable to better understand how CNTs with different helicities nucleate and grow through 
their synthesis process, in the hope of controlling helicity. 

In a preceding paper, I proposed a simple model of "nucleation and growth of single wall carbon 
nanotubes" [T]. I showed that it is possible to build embryos of CNTs on a graphene plane just by 
adding a number of bi-interstitials, and to make them grow by a similar mechanism, involving a very 
restricted number of defects, namely carbon interstitials. Stone- Wales defects and dislocations. 

In the present paper, I continue this work, which was done systematically for a lot of geometries 
and helicities, and I show that in the framework of the model some helicities are clearly favored against 
others. I compare these results with experimental data, which are discussed. 

2 Model and calculations 

The model for nucleation and growth of CNTs used here corresponds to the so-called root-growth 
mechanism. A good example of the use of this mechanism is the vision proposed by Maiti et al. |2]: 
a graphene layer wraps metallic catalytic particles oversaturated in carbon. Embryos of CNTs, fed 
with interstitial C atoms dissolved in the catalyst, can form on this layer, and subsequently grow from 
their foot. 

More recent papers go in the same direction. For instance, Gavillet et al. [3J suggest a root-growth 
mechanism, after a study both by high-resolution transmission electron microscopy and quantum- 
molecular-dynamics simulations. They state that they "will show that carbon atoms can be added at 
the root of a growing tube by a diffusion-segregation process occurring at the surface of the catalytic 
particle." 

Similarly, Ding et al. in the framework of the VLS model (vapor-liquid-solid), describe the 
tube growth "by small graphitic islands that form on the supersaturated cluster surface." 

The growth model suggested by Nasibulin et al. [5J is again quite similar. A carbon layer is said 
to be present on the surface of the metallic catalyst particle. In their fig. 5, they suggest that the C 
interstitials can come both from the saturated catalytic grain and from carbon entering continuously 
the grain through the opposite face of the particle. 



Figure 1: Nucleation of an embryo of (12,0) CNT on a graphene plane: 24 C interstitials are added 
in the 12 cells delimited by the dark lines and atoms. 

On this basis, my model lies on the fact that an embryo is nucleated in the graphene plane, so 
that it can grow from the foot, perpendicularly to the plane. Negative curvature on the graphene 
sheet comes from the presence of heptagonal rings, and positive curvature from pentagonal rings. The 
cap of the embryo must contain exactly 6 pentagons, with exactly 6 heptagons around the foot of the 
embryo [1]. 

To build such an embryo, let us consider, on the honeycomb network of a graphene plane, an 
hexagonal (regular or not) contour of elementary hexagonal cells, containing such hexagonal cells, 

being an integer greater or equal to 6. We add 2N C interstitials, two per cell, coming from the 
carbon-saturated metallic catalyst. An example of this process is given in fig. (T) a 12-cell contour is 
delimited by the two dark lines, 24 C interstitials are added, and bounds are rebuild in such a way 
to get a (12,0) zigzag tube. There is a lot of possibilities to build bounds between the interstitials 
as well as between them and the other C atoms of the plane. Among all these possibilities, there 
are several arrangements giving rise to an embryo with the required 6 pentagons and 6 heptagons. 
Such arrangements can be labeled by the usual indices n and m, corresponding to a (n, m) tube with 
n + m = N . One can see quite easily that by such a process all the helicities, i.e. all the n and m 
values obeying n + m = N , are reachable; moreover, in general, there is a number of isomers for each 
helicity. We can recall here the well-known fact that all helicities can be described by taking both m 
and n as positive or zero integers, with the condition m < n (which implies m < N/2). 

Once the interstitials are introduced and the bounds established, it is necessary to let the bond 
length and angles relax towards equilibrium and to check the stability and energy of the proposed 
defect geometries. I performed energy minimizations by using the Tersoff model [6]. The embryos are 
built on a parallelogram-shaped portion of a graphene plane, with periodic conditions in the x and y 
directions. This gives rise in fact to an embryo hexagonal network; I made calculations with varying 
the distance between embryos, in terms of carbon hexagonal cell rows, from one to three rows. 

3 Results 

I made calculations for embryos with N = n + m equal to 12, 15, 18 and 20, and for all possible 
helicities (from zigzag to armchair) in each case. In the following, all energy values are given relatively 
to perfect graphene with the same number of atoms. On fig. [21 I give the energy values for these 
geometries versus the m index, which is such that m = for zigzag tubes and m = n for armchair 
ones. Both total energy and energy per added atom values are given. 
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Figure 2: Energy of embryos of CNT versus m. N = n + m is 12 for white dots, 15 for diamonds, 18 
for triangles and 20 for black dots, (a): total energy; (b): energy per added atom. 
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Figure 3: Energy of embryos of CNT versus chiral angle. 

On fig. [3l the same energy values are given versus the chiral angle, whose value 6 comes from the 
relation 6 = arctan(-v/3 x m/[m + 2n)). 

On fig. m the same energy values are given versus the theoretical tube diameter, whose value D 
comes from the relation D = aV + + nm j vr with the lattice parameter a = 0.249 nm. 

For illustrating the geometry of the obtained embryos, I show on fig. [5] an example of tubes with 
= 12 for three helicities. 

I mentioned above, in section [21 that there is in general a substantial number of isomers for each 
helicity. This is especially true for values of n + m = larger than 12, when there are several possible 
configurations for the 6 heptagons in the hexagonal contour of the embryo. In practice, there is a 
small distribution of energies for the different isomers of a given (n, m) embryo, of the order of 1 eV. 
One notes that for zigzag embryos, there is only one isomer, when the 6 heptagons are situated at 
the vertices of the hexagonal contour; this is true for all sizes of zigzag embryos, and does not depend 
on the fact that the contour hexagon is regular or not. For all other helicity cases, different isomers 
exist in general, corresponding to different possible repartitions of heptagons on the vertices or sides 
of the contour. An example is shown in fig. [6] for an embryo of a (9,9) armchair nanotube, with three 
possible configurations of the heptagons on the sides. 
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Figure 4: Energy of embryos of CNT versus tube diameter. 




Figure 5: Tube embryos with N = 12 for three hehcities, from top to bottom: (6,6), (9,3) and (12,0). 
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Figure 6: Nucleation of three isomers of an embryo of (9,9) armchair CNT on a graphenc plane: 36 C 
interstitials are added in the 18 cells delimited by the external lines. The 6 heptagons are situated at 
different places on the sides of the hexagonal contour. 
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Figure 7: Growth of = 12 CNTs. black dots: initial embryo; triangles: one row added; diamonds: 
two rows added; white dots: three rows added. 
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Starting from the calculated embryos, it is easy to grow the CNTs, row by row. For each row, it 
is necessary to add between the foot and the cap 2N carbon bi-interstitials. I did such calculations, 
up to three rows, and got the energy values corresponding to the obtained geometry after the Tersoff 
minimization. An example is given fig. [71 

4 Discussion 

It seems very obvious from fig. [2] and [3] that the energy of embryos is always maximum for m = 
(zigzag tubes) and is minimum for values of m close to n, that is tubes close to the armchair geometry. 
The energy difference between the extreme geometries is about 2 to 4 eV, corresponding roughly to 
0.1 eV per C interstitial. The fact that there are several isomers in the armchair geometry whereas 
there is only one for zigzag tubes goes in the same direction: armchair or quasi-armchair tubes are 
predicted to be favorised in such a nucleation and growth process. 

There are several papers in literature dealing with CNT production and chirality measurement. I 
review now some of these papers, showing that the conclusions are quite far from unanimity. Cowley 
et al. [7] observed in CNT ropes produced by laser irradiation that the predominant helicity was 
armchair - (10,10) - or close to it. Qin et al. [8j found, also in laser-produced tubes, that in most 
cases the helical angles were evenly distributed, except a few cases where armchair structure was 
preferred. Bernaerts et al. [9] worked also on laser tubes and concluded that the ropes comprised 
largely tubes of armchair chirality, but not exclusively (10,10). Henrard et al. [TU] found that the 
situation was somewhat more complex; pointing the disagreement between the preceding papers, they 
claimed that within a given bundle no chirality was favoured. Colomer et al. |1H I12j worked on small 
bundles produced by catalytic chemical vapor deposition (CVD) and claimed that except for a weak 
tendency for the "armchair structure" , no preferred helicity was found. Bachilo et al. [13] observed a 
lot of semiconducting tubes grown in high-pressure carbon monoxide (HiPco) and obtained a random 
distribution of helicities; in a subsequent paper pi], they analyzed also tubes obtained by CVD and 
found then a significantly stronger preference for near-armchair structures. Li et al. [151 [T6] published 
a strong preferential growth of semiconducting tubes in the case of a plasma enhanced CVD technique, 
very different than in the HiPco or laser processes. Gavillet et al. [T7] concluded that "the helicities 
of the different tubes in a bundle obtained from the arc and the laser methods are distributed at 
random" and suggested that in the case of CVD, some bundles contain tubes with identical helicities. 
Vitali et al. [18], in the case of HiPco tubes, did not find any preferential helicity in the tubes that 
they analyzed. Miyauchi et al. [19] studied tubes synthesized from the alcohol catalytic CVD method 
and claimed that chiralities were distributed predominantly close to the armchair structure. Meyer 
et al. [20], working on CVD tubes, also found quite dispersed helicities. It is worth mentioning the 
extended work of Liu et al. [21| who studied 124 single-walled CNTs synthesized by the arc method, 
and noticed a slight preference for the helicity 15°-30°, quite close to the armchair structure. 

Reich et al. j22l ^3] gave a comprehensive theoretical study of the nucleation and growth of CNT 
caps on a catalyst surface. They pointed the fact that in a root-growth mechanism, the initial geometry 
of the cap controls the tube chirality. In such a framework, the epitaxial relation between the cap and 
the catalyst surface is of major importance for controlling chirality. 

An interesting work is the one by Lolli et al. [23]. They proposed a way to vary the distribution 
of (n,m) tubes with high specificity, by using CoMo catalysts and two C-rich gases, CO and CII4. 
With CO feed, they got rather small tube diameters, with a strong preference to the armchair or 
quasi-armchair geometry, like Miyauchi et al. With CII4, by contrast, tubes appeared near both the 
zigzag and armchair lines. 

In a recently published book. Page et al. [25j gave an extensive review on the mechanisms of 
single-walled CNT nucleation, growth and chirality-control, mainly focused on quantum-mechanical 
molecular dynamics simulations. In the experimental part of their review, they described not only 
the classical techniques like arc discharge or CVD but also more non-traditional methods using new 
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catalysts such as Si02, SiC and AI2O3. 

Finally, Sankaran |26] presented the recent advances and potential of plasma technology for CNT 
chirality control. He pointed that the catalyst structure can influence the nucleation of specific chi- 
ralities. Varying the plasma process temperature could also be a way to controlling the chirality. 

It is hard to deduce a precise tendency from this literature: for the different CNT growing tech- 
niques, is there a preferential nucleation of some helicities? The cited papers are somewhat contra- 
dictory on that question. However, it seems to be an overall preference for armchair-type CNTs; see 
for instance [TJ |9l [T^ [19] and particularly |24j . This tendency to get CNTs more on the armchair 
structure side is coherent with the findings of the calculations exposed in the present paper. 

I wish to add a few words about energy calculations. Li et al. [161 fig 4a] show density functional 
theory calculations of the formation energy of various CNTs versus tube diameter d: for tubes without 
foot or cap, they get a variation; for d = I nm, E is roughly 0.1 eV/atom. They also notice a 
small preference for semiconducting tubes. Liu et al. [21J, using ab initio calculations, found less than 
20 meV/at of energy difference between the different CNT structures. Reich et al. [22. table I] gave 
a lot of ab initio data on (10,0)-tube energies: 0.137 eV/at for the tube without cap, and between 
0.29 and 0.37 eV/at for the tube cap. My own calculations give 0.08 eV/at for a (12,0)-tube growth 
(fig- El) and 0.43 eV/at for the same tube cap nucleation. 

5 Conclusion 

The present paper develops the calculations for nucleation and growth of CNTs from a graphene sheet 
at the surface of a catalyst. The geometrical constraints of this process induce some energy differences 
between the different tube geometries, which are clearly in favor of armchair or quasi-armchair tube 
structures. 

There is no totally unambiguous experimental evidence of such a preference in the classical tube 
growth methods. However, several studies notice such a tendency to nucleate preferentially tubes of 
armchair-type chirality. 

References 

[1] F. Beuneu, Sol. State Comm. 136 (2005) 462. 

[2] A. Haiti, C.J. Brabec, J. Bernholc, Phys. Rev. B 55 (1997) R6097. 

[3] J. Gavillet, A. Loiseau, C. Journet, F. Willaime, F. Ducastelle, J.-C. Charlier, Phys. Rev. Lett. 
87 (2001) 275504. 

[4] F. Ding, K. Bolton, A. Rosen, J. Phys. Chem. B 108 (2004) 17369. 

[5] A.G. Nasibulin, P.V. Pikhitsa, H. Jiang, E.L Kauppinen, Carbon 43 (2005) 2251. 

[6] J. Tersoff, Phys. Rev. Lett. 61 (1988) 2879. 

[7] J.M. Cowley, P. Nikolaev, A. Thess, R.E. Smalley, Chem. Phys. Lett. 265 (1997) 379. 

[8] L.-C. Qin, S. lijima, H. Kataura, Y. Maniwa, S. Suzuki, Y. Achiba, Chem. Phys. Lett. 268 (1997) 
101. 

[9] D. Bernaerts, A. Zettl, N.G. Chopra, A. Thess, R.E. Smalley, Sol. State Comm. 105 (1998) 145. 
[10] L. Henrard, A. Loiseau, C. Journet, P. Bernier, EPJ B 13 (2000) 661. 
[11] J.-F. Colomer, L. Henrard, Ph. Lambin, G. Van Tendeloo, Phys Rev B 64 (2001) 125425. 



7 



[12] J.-F. Colomer, L. Henrard, Ph. Lambin, G. Van Tendeloo, EPJ B 27 (2002) 111. 

[13] S.M. Bachilo, M.S. Strano, C. Kittrell, R.H. Hauge, R.E. Smalley, R.B. Weisman, Science 298 
(2002) 2361. 

[14] S.M. Bachilo, L. Balzano, J.E. Herrera, F. Pompeo, D.E. Resasco, R.B. Weisman, J. Am. Chem. 
Soc 125 (2003) 11186. 

[15] Y. Li, D. Mann, M. Rolandi, W. Kim, A. Ural, S. Hung, A. Javey, J. Cao, D. Wang, E. Yenilmez, 
Q. Wang, J.F. Gibbons, Y. Nishi, H. Dai, Nano Lett. 4 (2004) 317. 

[16] Y. Li, S. Peng, D. Mann, J. Cao, R. Tu, K.J. Cho, H. Dai, J Phys Chem B Lett 109 (2005) 6968. 

[17] J. Gavillet, J. Thibault, O. Stephan, H. Amara, A. Loiseau, Ch. Bichara, J.-P. Gaspard, F. 
Ducastelle, J Nano Nano 4 (2004) 346. 

[18] L. Vitah, M. Burghard, M.A. Schneider, L. Liu, S.Y. Wu, C.S. Jayanthi, K. Kern, Phys. Rev. 
Lett. 93 (2004) 136103. 

[19] Y. Miyauchi, S. Chiashi, Y. Murakami, Y. Hayashida, S. Maruyama, Chem. Phys. Lett. 387 
(2004) 198. 

[20] J.C. Meyer, M. Paillet, T. Michel, A. Moreac, A. Neumann, G.S. Duesberg, S. Roth, J.-L. Sauva- 
jol, Phys. Rev. Lett. 95 (2005) 217401. 

[21] Z. Liu, Q. Zhang, L.-C. Qin, Phys Rev B 71 (2005) 245413. 

[22] S. Reich, L. Li, J. Robertson, Phys Rev B 72 (2005) 165423. 
[23] S. Reich, L. Li, J. Robertson, Chem. Phys. Lett. 421 (2006) 469. 

[24] G. Lolli, L. Zhang, L. Balzano, N. Sakulchaicharoen, Y. Tan, D.E. Resasco, J. Phys. Chem. B 
110 (2006) 2108. 

[25] A.J. Page, Y. Wang, K.R.S. Chandrakumar, S. Irle, K. Morokuma, in Electronic Properties of 
Carbon Nanotubes, edited by J.M. Marulanda, InTech (2011) p 521. 

[26] R. M. Sankaran, J. Phys. D: Appl. Phys. 44 (2011) 174005. 



8 



